Soliton solutions of the classical lattice sine-Gordon system 



B. Enriquez 

Abstract. We study the soliton-type solutions of the system introduced by B. 
Feigin and the author in [EF]. We show that it reduces to a top-like system, and we 
study the behaviour of the solutions at the lattice infinity. We compute the scattering 
of the solitons and study some periodic solutions of the system. 



Introduction. 

In this work, we study solutions to the lattice sine-Gordon system introduced by 
| B.L. Feigin and the author in [EF]. This system consists of two families of compatible 
flows, analogues of the sine-Gordon and mKdV flows. We recall that the solutions 
,-0 ■ studied in the continuous context are characterized by the requirement that they be 
. stationary with respect to a linear combination of those flows. We study the analogous 
^ ! problem here, and reduce it to a top-like system. We then show how to pass from 
the group variables satisfying this system to the local variables (inverse scattering) ; 
in particular, we find solutions with variables tending to a constant at lattice infinity 
(Sj ■ (soliton solutions). We then study the scattering of these solutions, and we find a 

purely elastic behaviour. Two features are rather different from the continuous case : 
the absence of charge of the solitons and their non-fermionic character (solutions with 
multiple poles exist). This is due to the difference between the relevant subalgebras 
of sl<2,: positive Cartan and negative principal subalgebras instead of positive and 
*7h 1 negative principal subalgebras in the continuous case. 

Finally, we study some periodic solutions of the system. We find that such 
solutions are characterized by the condition that the difference of the two points at 
r-| ! infinity on the spectral curve is torsion in its Jacobian. 

We express our thanks to B. Feigin for discussions related to this problem and 
to Mme Harmide for typing this text. 

1. Review of continuous sine-Gordon system. 

Let 4> be a function on R and let us consider the matrix M X (X) = Pexp (A+ 



(j)h).n x , with n x = n x exp ^ J ±0 ° UiA l , and d x +A+(j)h = n x (d+A+ w iA l ) n x 1 i 

i>0 i>0 
odd odd 

/i=^J , A = ^ J^j ([DS]) ; Pexp J^ qo (A + <ph) is developed around A = 

and takes the form 

/ 1 + A e 2 ^) e~ 2 - + ■■■ Jle^... W 

V AP e~^ + --- 1 + AP e~ 2 ^ f y e 2 ^ + ■ ■ ■ V 

v J— oo J —oo J —oo ' v 

Then the mKdV flows are d^ dV M x = M x A 2n+1 , n > and the sine-Gordon flow is 
d$ G M x = A _1 M X . We introduce higher sine-Gordon flows by d^ G M x = A~ ( - 2n+1 ^M x . 
We have 

d x dl G ^(x)=e 2 ^ -e~ 2 ^\ 
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d x df G <p(x) = e 2 ^ [ X e~ 2 ^ f (e 2 * - e - 2 ^)- e - 2 ^) f e 2 ^ f (e^-e" 2 ^), 

J —oo J —oo J — oo J — oo 

d x d s n G ^(x)=e 2 ^) f e- 2 ^d y d s n G Mv)-e- 2 ^ x) f e 2 ^d y d^Mv)- 

J — oo J — oo 

These equations are well defined for functions satisfying the boundary conditions 

e 4<p(x) — i faster than any inverse polynomial, and they preserve these bound- 
cc^— oo 

ary conditions. If (p is such that in addition e 4 ^) — ► 1 faster than any inverse 

x^-\-oo 

polynomial, and if the conditions 




are met (they are a rewriting of [A, Pexp (A + <f>h)] = 0), then the higher 
SG equations can be rewritten with +00 replacing —00, at point ip ; these con- 
ditions are preserved by the flows. [Note that any solution to this system, such 
that e 4lfi — > 1 and d^ G p — > 0, faster than any -^j- when x — > ±00 will also sat- 
isfy the conditions f_ (e 2ip — e~ 2(p ) = • • • = 0, since we deduce from d^ G (ph) = 
(Ad- 1 (Pexp/" oo (A + ( /)/i))(A-( 2fc + 1 )))^ A o that Ad _1 (Pexp /^(A + ^/i))(A _1 ) has 
diagonal elements equal to zero, and from <9f G (l + </>A -1 ) = degree —1 part of 
Ad- 1 (Pexp/^ oo (A+ ( /)/i))(A-( 2fc + 1 )) we deduce that Ad _1 (Pexp /^(A+^XA" 1 ) 
is proportional to A -1 ; since its square equals A -2 , it is equal to A -1 . These 
two equations are consequences of d^ G M x = M X (M~ X A~( 2k+1 ^ M x )> and d^ G n x = 
(M" 1 A-( 2fc + 1 ) M x ) <0 n x . The first example of this statement is /^ (e 2 ^-e" 2 ^) = 0; 
it is a consequence of d x dQ G p> = 8q G [p(oo) — p( — 00)] = (conservation of topo- 
logical charge).] 

The higher sine-Gordon flows commute with each other and with the mKdV 
flows. The solution of this system is formally written 

Mx(X; t s n G , t« dv ) = exp(£ 4 G A-( 2 ^) 9o exp( £ fJ^A 2 "* 1 ). 

n>0 m>0 

Then to extract e v from this matrix, we have to take matrix coefficients of this group 
element in some integrable representation. It means that M x {\) has to be replaced 
by a variable M X (X) living in the central extension of the loop group, obeying the 
same equations as M X (X). Now the Cartan part of M X (X) is exp(tph) times a central 
term : it will be equal to (A |M x (tf G , t ] ^ l dV ) |A ) if |A ) is a level-one highest weight 
vector, such that h\A ) = 0. Now, if |Ai) is a level-one highest weight vector, with 
h\Ai) = I Ai) , we have 

e *<*5°,*£* v ) = (Ailg^tg ,^)^!) 

(Ao|M x (^^W)| Ao) - 
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(cf. [OTU]). 

In the case of the soliton solutions, we take go = e^ lF(yZl ^ ■ ■ ■ e ( ^ nF ^ Zn \ with 
[A»F(z)] = z n F{z). Let n(t s n G ,t« dV ) = (A 1 |M a; (t^ G , t£ l dV )\A 1 ), then d s n G r % = 
e^+ElLi J§S because [A"» <7o] = £" =1 ^ SO ' ^ G (^) = £? =1 sp^M ; 
similarly, d n (e^) = £™ =1 ^ggrCe*')- This shows : 

Proposition. — The mKdV and higher SG flows span a finite dimensional family 
of vector fields on the multisoliton solutions. 

Note also that go can be considered as a monodromy matrix ; from its expression 
nlLi ex P(£fcez hh k z~ ) we deduce that as a function of A it has poles at A = Zi. 

2. The lattice sine-Gordon system. 

Let as in [EF] x^ y^i G Z) be a system of lattice variables, with Poisson brackets 
{xi, xj} = XiXj, {yi, yj} = y^yj, i < j, {x h yj} = -Xiyj, i < j. Consider the matrix 



g = n 



' 1 OUl Vi \ f 1 0\ 1 1 [»_,+ 



%= — oo 



Az, 1/VO 1 7 ' V l/Xx + l/y 1 --- 1/10 1 



n Ax t Vi n. 
*>1 Axi + l/j/i- i/i+l/Aa;i + - u 

n (TT Ax t ^ \-i 

u UIi>l Axi+l/j/i- j/i+l/Axi+.-J 

The lattice sine-Gordon and mKdV flows can be expressed as 

(1) af G 9 = A- <2 " +1 > S and dZ* v g = gh\",n>0 ; 

hereA=(° J ). Let M be the matrix ^ J) (J 

Proposition. — The following relation is satisfied 

{M®,g} = p(M®g) , 

(0 A* A 
J ' ^ = 

0\ , _ /A* 0\ 
A* J ' fti ~ V " A V 

Proof. For the matrix elements of M corresponding to the simple roots 
(E-oo^ and E-oo^)> this statement is {E-oo x ;>#} = - 2e -i#> {E-oo^'^l = 
-2/o0, and it is a consequence of {E-oo^^}' = - 2 (fi ,e -ifi'~ 1 )-fi , 5 {E-oo^'^}' = 
—2(gfog~ 1 )~g which follow from [EF] (here {a, 6}' = {a, 6} — (dega)(deg&)a&, degXi 
= 1, degyi = —1). We then check its compatibility with {M®,M} = r L (M <g> 
M) — (M ® M)r R : it is (5 ® l)(p) = [p 23 ,p 13 ], (where = [r, x <g> 1 + 1 <S> x]). 
Recall that p G b + <g> &_. Let £,77 G 6_, this equality means that (p, [£,77] ® 1] = 



[(p 13 , (p 23 , rj 2 )] ; since [p, £<g)l] = £ this equality is true (we have used the duality 
between b + and &_)• ■ 
Let us set Hi = constant coefficient of tr(A~( 2 * +1 )M). The family Hi is Poisson 
commutative. Then 

{Hug} = J2- 2d i-j(f-j9) -2a l . J (e I l J _ 1 g) - (c^j - , 
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where we have set M = ^ k(A) ^{X) ) an< ^ a ^ = ^»>o ai/ ^' etc --- 

The Hamiltonian flows generated by the Hi commute with the mKdV flows 
(since Hi are in the Poisson algebra generated by b\ and Co, which commute to the 
integrals of motions) ; these flows will be denoted <9#\ . Generally they will make sense 
only at the level of formal variables (whereas d% G = fb n + e^ n _ 1 will have concrete 
realizations), because a;_j and di-j will be infinite. Nevertheless, for solutions of the 

system such that x n , y n — ► a, we can consider as we did in the continuous case that 

±00 

cii — di = bi — Ci-i = and that the linear spans of dn 1 , • • • , <9h„ and df G , • • • , d% G 
are the same. 

Let us now show how the vector fields d^ G operate on the space of lattice vari- 
ables, rapidly tending to a at infinity. Let us express those fields as 

df G {\nx n ) = -2^2yi+2^2xi+x n = (x n -y n _i)-(y n _ 1 -x„_i) + (x n _i-y n _ 2 ) • • • , 

i<n i<n 

df G (lm/ n ) = (y n - x n ) - (x n - y n -i) H , 

and (lnx n ), d% (hay n ) are given by (1), where the matrices 



Ad- n 



1 0\ (1 -yA fO A" 1 



-\xj 1 j \ 1 j \l 



1 -yA (I (A/0 A" 1 
1 J \ -Xx j+1 1 J \ 1 



Ad- n 

j=-oo 

are understood as 

i 

(3) lim Ad -1 TT 

j=-N 

and 

(4) 

1 °\ tj f 1 -vA ( 1 A 1 (-, a- 1 



1 0\ A -yA 1 (% A" 1 

-Ax, 1; \0 1 J f^KA -f 



iv^oo^A-A^ iJ.n^U i J ){-\ Xj+1 i)'r^\i 
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(these limits are well defined, because of the assumptions on x u yi ). Thend££(l na >)) 

is a linear combination of J^n„,^o4 , h ^ 4 _l„ xf 1 ^ • • • x- £2fe ' ) 9 1 SG lnx- e2fe+1 ' ) 

where = x n for e = 0, and y n for £ = 1, and fcj — ki-i = or l(ko = 0). 

These equations are well defined for sequences satisfying x n , y n — > a faster 

n^ — oo 

then any and they preserve this boundary condition. If (x n ,y n ) satisfies the 
same conditions for n — > +oo, and if all coefficients of A fc in An,n' +o,Bn,n' — Dn,n' 
and -B^ at' — ACjv,at', tend to zero faster than any jjr-^r when N,N' — > +oo, 
Ann' C nn > \ _ ( 1 \ / 1 y N > 
B N>N t D NN > J \Xx- N 1J \0 1 
of equations can be rewritten with +oo replacing — oo ; these conditions on the lattice 
variables are preserved by the flow. Once again, any solution to the system, such 
that x n , y n — > a and d tk x n , d tk y n — > 0, faster than any for x — > ±oo, has to satisfy 
the above conditions on A N N > etc... Indeed, the system has to take the forms 



where M/v jv' = 



then the system 



-f A -1 



n 



for n = 



1 l/{\Xi 

1 



i • • • 





- 1 1 



o r 1 )' a,,<1 



§ A" 1 



for n' = 



1 

i/yi 



This implies that (Mjv,jv) 

1 a\/i C 



aAC -> 0, and ( Q 1 g £> M I 



or 2A+aB — > 0. So we have M^ 1 ^, | , - ,, | .: / v . \ 



— — \ — * 

7 a i 



any inverse polynomial ; the conditions (2) are satisfied. 



OA (l —[Xxi-i ■ ■ \ ft' 

i ) \o i ; V° 1 

if-i x ~'\ M ( A c 

1 f J M N,N' = V 

1 -a\ M' C" 



for |i| large enough, 
should satisfy 2A + 



, should satisfy 2 A! = aB' , 



— — A - * 

a ) fast cm- than 
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3. Finite dimensional orbits of lattice flows 

Let G be the completion of S£2(C((X))), equal to the product 7r _1 (.B_|_) x 
Tr-^iV.),^ : S£ 2 (C[[X]]) -> S£ 2 (C),X h-> and 7T : ^(C^A" 1 ]]) -> S^C^A" 1 i-> 
0, and let A + be the subgroup of 7r _1 (i?_|_) corresponding to the Lie algebra a + = 
©fc>oCA fc (e+A/), and ii_ the subgroup of 7f _1 (iV_), corresponding to h- = Q) k >oChX~ 

There is a bijective correspondence between elements of A + \G/H_ and systems 
of variables Xi, yi, M^,i G Z, xi, yi scalars and elements of Lie 7r _1 (£?_)_), conju- 
gate to A _1 e+/ and such that M+ = Ad ^ J ^jjfr, ^i+i = Ad ( ^ l ) M ^ ' 

To see it, we express A + gH_ as A + n + ri-H_ and write n_ = 1/A^o + ^ 

• ^ ^ ^ ^j, and A + gw ± (k)H^ = A_|_n^(/c)n^(/c)i/_, and then pose 

= Adn|(/c) _1 (A _1 e + /) [where w + (k) = (w wi) k , w~(k) = (w wi) k w ]. It is 
clear what the images of the natural flows on A \ G/H_ are by this correspondence. 
Finite dimensional orbits of these flows are also in correspondence ; they are also the 
orbits where a linear dependance condition between some of these flows is satisfied. 

Solutions (x n ,y n ) to (1), (2), satisfying the conditions x n (t),y n (t) — > a faster 

n^oo 

than any for any fixed t, yield solutions to the above described system on 
(xi,yi, M^) : it suffices to identify Mf and M~ with the matrices defined in (3) 
and (4). 

We will now describe the finite dimensional orbits of A + \ G/H-, and describe 
those satisfying conditions at infinity (the soliton-like orbits) . 

Let g G G. It will lie in a finite dimensional orbit, iff a condition of the 
form (A" 1 + Efc>o a fc A2fc+1 )f = 9j2 k < N b khX k (b N ^ 0) is satisfied. Writing 
g = n+b-, this means the existence of some x G s^2(C[A, A -1 ]), such that (A -1 + 
^2k>o ttfcA 2fc+1 )n-|- = n+x, and xb- = b-^2 k<N b k hX k ; x has to be of the form 
A _1 + elements of homogeneous degree > and < 2N + bNhX N . The vector fields 
(1), (2) induce on x the flows 
(5) " 

xX n 



d s n G x 



X~ n 



and d* dV x 



bNX N + J2 k<N hX k 



,x 



n > 1, keeping a k and b k fixed. 

Let us indicate how to obtain (xi, yi, M^) if we know x(t). Suppose that g(t) cor- 
responds to the dot Xi (resp. yi). ThenAd^J ^ x(t), resp. Ad ( ^ ^\ x 

gives the matrix x(t) corresponding to the previous dot, yi-i and Xi-\ respectively. 
That these matrices again take the form allowed for matrices x(t) determines yi-\ 
and Xi-i : posing x(t) = ±b]yhX N + a^-ieA^ -1 + cwfX N + lower degree terms at 
dot Xi (resp. yi), we deduce 

Vi-i{t) = - — -(£), resp. Xi-i(t) = — — — (£) . 
cat ajv-i 
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The values of the other variables Xf-(t), Vk(t) can be defined inductively. can now 
be obtained as follows : Mf = n^ 1 A~ 1 n + , so M^(t) = - ^ — 1 - ^ fc+1 x(t), if g(t) 

corresponds to the dot Xi (resp. yi). 

In general, the system (5) can be solved noticing that x(t) follows an isospectral 
evolution, and that it corresponds to linear flows on the Jacobian of the curve det(p — 

x(X)) = fi 2 - (E k <NhX k ) 2 = V 2 - A- X (l + E k >o a ^ k+1 ) 2 = 0. If this curve is 

smooth, its Jacobian is compact and the solutions (xi, yi, M^) are expected to have 
a quasiperiodic behaviour. In any case, we will see that such solutions cannot have a 
soliton-like behaviour; such solutions will correspond to rational curves with double 
points. 

4. Soliton solutions of the lattice flow 

4.1 Integration of the system 

Let us determine now the (x(t),ak,bk) corresponding to a constant solution 
Xi = yi = a. We have for such a solution, N = so at dot Xi, x = A -1 + b\h and 
Xi(t) = — 2bi, and at dot yi,x = A -1 — b\h and yi(t) — —2b\ ; so b\ = — | ; the 
constant solution corresponds to the rational curve fj 2 = (|) 2 + A -1 . 

Let us suppose that (x^, yi, M^) has soliton behaviour, and corresponds to a fi- 
nite dimensional orbit. Then by (3), (4), we have for any fixed t = M^it) — > (1+ 

i — > — oo 

^) -1 / 2 y ^ _ s _ a J (this means that all coordinates of M^{t) in the natural coor- 
dinate system of Si2(C((X))) has the limit indicated). This implies that ^ — - — r^rr - 
' x (xi),( yi )(t) have the same limit ; X( Xi )(t) and X( y ^(t) denote the functions x(t) rela- 
tive to dots Xi and yi. We have now x^ Xi )(t) — > (J+^)2\y/2 (? a j ? an d since 

X( x .)(t) has no terms of homogeneous degree > N, we deduce that (j+f^yxyT^ * s 
a polynomial of degree AT, n^Li(l + 7»A). So the spectral curve has equation 

N 

(6) ^(A^ + ^IIa+^A) 2 ; 



it is a rational curve. This phenomenon is also observed in [NMPZ], II, 10,1. Pos- 
ing (A -1 + (f ) 2 ) rL=i(l + 7iA 2 ) = A -1 + Xli=o ^A 1 , we get the linear dependance 
condition df G + £ 2 =i did? dV = on our solution. If we assume all ct^'s to be real, 
the checking of Xi(t) — > a, yi(t) — > a for ti — > ±oo, other t^s fixed will be reduced to 
a problem of a finite number of flows. (Some solutions with complex a^'s have the 
same property, e.g. 1-soliton with Arg 7 e irrational.) From now on we suppose all 
cti's real. 
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Let us solve (5) in this case. Let x(X) = ^ ^A) ) ' an< ^ = ni^=i (1 — 

jM ; then a(A^) 2 = fi(Xi) 2 , so a(A^) = /Uj, (Aj, being some point of (6), lying over 
Aj. We deduce 

N N AT AT 

(7) a(A) = e| J] 7* II( A " ^ + E * II I^TT' £ = ±X • 

i=l i=l i=l j'=i 1 J 

Let us consider instead of the flows (5) their linear combinations 

d? k x=[(\- k x)+,x],k>0. 

Note that all df x are for n < N, and for n > N they will all be proportional to 
x]. The vector fields 9 t ' induce on 6(A) the flows ^d t >b = (X~ h b) + a — (A~ fc a) + 6 

(here (A*)+ = X l l i>0 ). Evaluation at Aj gives £ JX,^ 1 - = /M A 6 )+( A ») = 

h(x) k+l 
HiXi- coefficient of A fc in -^f- = / u i (-l) fc+1 Yl \ }.\. > so 

ic^all different 

fl A - o (-l) fc+ V^, ^ 

*'* *~ n — (x -x) 2-*> a ^'" a 3n-h-i- 



all different 



Consider now the forms u s = X s ^-,s = —1, • ■ • , N—2, and the functions F 8 (t') = 
J2iLi J^^ T ''^ l ^ T '' uj s (Pq a fixed point on the curve). They obey 



V \s+l 



all different 



by a Vandermonde system computation. So F s (t') = —2t' s+1 + est. But 

Ffn= f^ / a< ™ = v r ~ 2(r2 - ^ )2)N ~ s ~ 2 

iiL ((t^+A-^^niiCi+TiA) ^4 nf=i(r 2 +7,--(t) 2 ) 



with r 2 = (|) 2 + A 1 ; more precisely r = l^/Y\iLi(^ + 7i A )- 

Assuming all 7^'s different, the fraction in r is decomposed as 



_j2 (-^) N - S - 2 ( 1 1 



t n^(7i - 7i) vet ) 2 - - ^ r + A/ifF- 7 ^ 

for z = —1, • • • , AT — 2. This gives 

m h v<«>"--*iw-ft-T»> w<i>"-T» " i+1 +1 



C s+ i are constants. 

If we pose C(A) = A" 1 Jlilil 1 - and a i = K^i)/ U^ii 1 + lj v i)-> we nave 
similarly 



4.2. One-soliton solutions 

Consider the case N = 1. 71 = 7 is real. We find, with e = — 1 



1 + Ke^^W^i rjT l + K 'e 2t 'oVW^ 

(9) T l = J(-y- 1 —— ===== , £71 = d(-y -7 — /TTTT^ 

X, X' are (complex) constants of integration, related by 

g _ § - a/(IF^7 

For the solution to have the properties Xi(t' ) — > — - , /Ui(t ) — > — - for t — > ±00, 
we need a/(§) 2 — 7 to be real ; so the condition is 7 < (f ) 2 . We note that the values 
of 7~i, <7i are independant of the sign of a/(§ ) 2 — 7- We compute 

(l-ife^') 2 , (1 -iiV°) 2 

(10) Xi_i(*o) = a- —777— K2 t „. ,yi-i(t' ) = a- „ ^72 ^777 

(1 — i i C / e t o )(1 — ^e r o ) (1 — Ke l o )(1 — ^e £ o ) 

with = 2^/(1)2 -7t . We have (*{,), 2/i_i(*o) -»■ a fo r *d ±0 °- 

Replacing £ by 1 in (7) amounts to replacing a by —a. (10) can be considered 

as a one-soliton solution to the system. 

The values of Xi(t' ), yi(t' ) are given by formulae (10), with K and K' multiplied 

by (7^) 2 : it means that the lattice translation of the solution corresponds to a time 

translation of 



( n ) A *o = m 



Note that up to time shift, complex solitons depend on one variable (Arg K) and 
real solitons on no variable. 
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4.2 Breather solutions. 

Let us consider now the case N = 2, 72 = 71 not real. Rather than solve (7) 
explicitly, let us make the following qualitative remarks. Complex solutions depend 
(up to real shifts in t' and t[) on two real variables (ImCo up to integral multiples 
of Re — ; N „ 27I " 7l / r and ImCi up to integral multiples of Re — , sn 27r , -). 

2V(f) 2 -7i(7l-72) F & F 2^(1)2-71(71-72) ' 

The conditions for the solution to be real are r 2 = fi, or T\ and r 2 real, so the 

Tj = y"^_i In— — ^ 2 \^ 7l : should be conjugate for i = 1, 2 and this means that Co 

and C\ should be real : again the real solutions depend on no additional variable 
(up to time shifts). Let us again consider a complex solution and let us study its 
evolution for large times with at' + f3t[ fixed (a, (3 real not both zero). Since T\ = 
V(f) 2 -7i(*o + 72*i), T 2 = 2y/(f ) 2 - 72(^0 + 71^), Re Ti and Re T 2 tend to 
infinity except if (1 — 72^)/ a/(§ ) 2 _ 72 is pure imaginary (this fixes one value of ^). 
Using 

(12) r 1 r 2 (l-e T 0-^/(^) 2 -7 l (l + e T 0(ri+r 2 ) + ((|) 2 -7 l )(l-e T 0=0,z = l,2, 

we see that ri and r 2 tend to d=\/ (|) 2 — 71^2 so 2^ and j/i tend to a. Let us study 
now the behaviour of these solutions w.r.t. lattice periodicity. The lattice translation 
corresponds to some (a priori complex) shift in times (t ?^i) l— * (^o + Ao,^ + Ai). 
Ao and Ai should be analytic functions in 71 and 7 2 ; in particular, let us compute 
them in the case where 71 and 7 2 are real (71 7^ 72). In a region 

|*o + 72*il < l*o + 7i*il , 

the solution Xi(t' ,t[) tends to x\ (t' Q + 72^1), where x\ (t' Q ) is the soliton solution 
with parameter 71 described before. In this situation, the lattice translation corre- 

sponds to the shift (11), so A +7 2 A X = , 2 In- — v , 2 1 ; we have the same 

V(f) Ti f+V(f) T2 

equality exchanging 71 and 72, so 
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A = ) ^ 



27i(-l)* J-TdM 



7i-72^ v/(f) 2 - 7. f + V / (f) 2 -7 l ' 

a _ 1 y 2 (~ 1 ) i 1d !-a/(F^ 

1 7i-72^ v/(f) 2 - 7, f + V / (f) 2 -7, ' 

We see that in the case 7 2 =71, Aq and Ai are real. Note that the times direction 
(Aq,Ai) corresponds to a direction where the solution tends to a ; we have then 
Xi+N^o^ t[) = Xi(t' Q + iVAo, t[ + NAi) — > a when N — > 00. One can easily show that 
this convergence in exponential in N. So the solutions found here (with 72 = 71) 
have soliton-like behaviour (in the sense of 4.1) and can be thought of as analogues 
of the breathers (despite the fact that solitons are not charged). 
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4.4 iV-soliton solutions. 

Let us turn now to the general case : we have p real poles 71, • • • , 7 P , (7* < (f ) 2 ) 
and q pairs of complex poles 71, 71, • • • , 7q, 7q ; all pairwise different. We will show 
that the corresponding solution has soliton-like behaviour, and that it corresponds 
to the scattering of p solitons and q breathers. Now N = p + 2q. Let us discuss the 
effect of a large translation in times along some vector (a ? • • • ; &n-i)- (8) can be 
solved by 
(13) 

N — / N 

E' ^ir^ MW-^E £ ^.•••tv 1 (*j-,+c,-.)=t, 

i=i T i + VUi V ^ J=1 4 ^ 

all different 



As long as Re 2y/(%) 2 - 7; £V =1 ^ 7t ai ' ' ' 7i« J _ 1 are all non zero, all 

all different 

Re(Tj) will tend to infinity, which means that Tj tend to a/ (f ) 2 — 7j (up to order) 
and x a (^))J/a(^) tends to a. 

Let us discuss now the effect of a lattice translation. It corresponds to some 
shift in times (tj) 1— > (^ + Aj), z = 0, • • • , N — 1. We again compute Aj by analytic 
continuation of the case where all 7j's are real. Fix some index i. In a region |Tj| <C 
|Tfc|, A; 7^ i, we have — > ±a/ (f ) 2 _ 7fc and rr Q (to, • ■ ■ , £jv-i) tends to x a (Ti + const). 
Proceeding as we did previously, we find 



N 

E 



2(-70 Ar - 1 - a ln f " 



In our case, A' a s are real. Since Re(2(-1)^(§) 2 - 7; Ejli E '«« fc ,y 7i ai -"7 



(a — /( i ) 2 — 7i \ 
ln -|-j-^== = 1, it is not zero, by the assumptions. We deduce again that 

for any (tj), x a (^) and y a (t'i) tend to a as a — > ±00 ; this convergence is exponential 
in a, as is easily shown. So the solution has soliton behaviour. 

Let us indicate the regions of large times where the solution does not tend to 
a : we have |Tj| <C \Tk\, |Tg + Tg|, for i < p, k < p and k^Q, £>p+l; and 
\T t +f e \ < |T fc |, |1> + for i < > p+l,£^ £'. They correspond to solitons 

and breathers with parameters 7^,2 = 1, • • • , N — 2. 

Let us compute now the phase shift of the particle 7$, due to the scattering with 
other particles. 

We find, when tj goes from a/(§) 2 - Tj to -a/(§ ) 2 - 7j, for j 7^ i , ATj = 

if i > p. So the phase shifts are the sums of two-particle contributions ; this is a 
situation of elastic scattering. 
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4.5. Multiple-poles solutions. 

Assume that the poles 7o,7j,7j(z < P, j > P + 1) have multiplicities a^ctj. 
Formula (8) is then replaced by 
(14) 



V(f) 2 -7^ n ^ + V / (f) 2 "7, 
t' s+ i + est, where AT = ^i=i because lim 7i ^ 7 Yh=i tt F f 7i) _ ^ 



f (7i) _ ^ (fc) (7) 

' a: ' 

which is solved by 



(15) {^-YY\n 3 1% =TjB , )9 = 0,...,a 4 _i, 

where the Tj j( g are linear combinations of the times t' s+1 . The effect of a lattice 
translation and the large time behaviour of this system are studied in the same way 
as for the simple poles solutions : so multiple poles again define solitonic solutions, 
whose scattering can be computed similarly. 

We note two main differences with the continuous case : absence of topological 
charge, and no fermionic nature of the solitons. 

5. Periodic solutions of the lattice flow. 

We now examine the possibility of finite dimensional orbits of the lattice flow 

yielding periodic (in i) solutions Xi(t),yi(t). Let v be this period. The condition for 

{ ^ * X"' 1 + ■ ■ 

periodicity is the existence of an element b in N + , of the form . , . v 

\ *A A -r ■ ■ ■ 

*' 7^ (for example), commuting to X( x .} (t). This condition is preserved by the flows : 
b has to follow the equations dffb = [6, ( 1+ ^ J Afc+1 x)+] = [( l+ ^ J Afc+1 x)-,b] and 

QtJ) = [( ^+£L;v^ fc ) - ,b] = WhiCh dearly Pr6SerVe 

its form. Writing b = P + Qx, P, Q rational fractions in A, the conditions on P and 

Q are : P polynomial, Qx polynomial and P 2 + Q 2 detx = 1. Thus the condition 
is the existence of a meromorphic function on the spectral curve, with only zeroes 
and poles at the points at infinity oo . The order of this zero (or pole) is v, and we 
see that we have z/(oo + — oo - ) = is the Jacobian of the curve; v is also the lattice 
period. Thus the periodicity condition is that the difference of the infinity points is 
torsion is the Jacobian of the spectral curve. 

We compute now the dimension of the space of ^periodic solutions among all 
spectral curves of degree N. Writing det x = —A -1 rii=i + (1 ~~ x~) an d Q = ^,R, S 
polynomials of minimal degree, we find that A should divide R. We find also that 
S should divide rii=i +1 (l ~~ x~)' an< ^ a ^ so x > but we can restrict ourselves without 
loss of generality, to the case where no factor of rii=i +1 (l ~~ X") divides x. Let now 
R — XRq, S = 1 ; the condition Qx polynomial is now satisfied, and we need only to 
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satisfy P 2 = 1+P^A n 2 =i +1 (l- The degrees of P and #0 are v and v-(N+l) ; let 
us write [An-5i +1 (l~ A)]' 72 [l + ^fiv+i) +0(\-^( N +^)} ; we obtain 

^ = (t A n?5i +1 (! - ir)] 1/2 )- + a^iv+ti » for p o the remainder of the division of P 
by P and P' = K{- W™? 1 A,)" 1 / 2 - The space of possible fractions has dimension 
2{v — N — 1), and 2u — (N + 2) equations must be satisfied in it ; this imposes 
iV conditions on the (Xi)i=i,---,2N+i- It would be interesting to see if some natural 
Poisson brackets on the (Aj) can be introduced, and what kind of submanifold the 
periodic solutions will be in the space of all solutions (or of hyperelliptic curves) with 
respect to this Poisson geometry. 
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